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Abstract 



Critical catalytic branching random walk on an integer lattice Z d is investigated for all 
d G N. The branching may occur at the origin only and the start point is arbitrary. The 
asymptotic behavior, as time grows to infinity is determined for the mean local particles 



5— i 

numbers. The same problem is solved for the probability of particles presence at a fixed 

lattice point. Moreover, the Yaglom type limit theorem is established for the local number 

of particles. Our analysis involves construction of an auxiliary Bellman-Harris branching 

process with six types of particles. The proofs employ the asymptotic properties of the 

(improper) c.d.f. of hitting times with taboo. The latter notion was recently introduced 

by the author for a non-branching random walk on Z . 
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^o '. 1 Introduction 

m 

• ■ Catalytic branching random walk (CBRW) on d- dimensional integer lattice is a model of par- 

ticles population evolution. We recall its main features. Each particle independently of others 
may perform random walk on Z d and produce offsprings at the source of branching located 
w.l.g. at the origin. Symmetric branching random walk (SBRW) on Z d studied earlier, e.g., in 
P], PJ and |2I1 is a particular case of CBRW on Z d (see [23]). 

The model under consideration was proposed in [20J for d = 1 and studied for other d E N 
in jl], [5], [13] and [23]. The analysis of CBRW in [13] and [23] has shown that similarly to 
many kinds of branching processes (see |12] ) CBRW on Z d is classified as supercritical, critical 
or subcritical. According to [23], the exponential growth (as time tends to infinity) of total 
number of particles in population and local numbers of particles as well is characteristic for 
the supercritical CBRW on Z d . The term local refers to the (number of) particles located at a 
lattice point. 

Quite different situation occurs for critical CBRW which is the main object of study in this 
paper. For example, for d — 1 or d = 2 the particles population degenerates with probability 1 
but survives with strictly positive probability for d > 3 (see [5], [7], [13] and [2D])- Moreover, 
the total number of particles conditioned on non-degeneracy has non-trivial discrete limit dis- 
tribution, different for d < 3 and d > 3 (see [I], [5] and [20]). Thus, in the model of critical 
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CBRW on Z d the asymptotic behavior (in time) of the total number of particles on the lattice 
depends on dimension d essentially and does not grow exponentially. As for the local particles 
numbers in the model of critical CBRW on Z d , earlier in [I] [2], [H], [13], [H] and [20] there 
were only established the asymptotic properties of the number of particles located at the source 
of branching. In particular, it turnes out that for all d G N the probability of the presence of 
particles at the source of branching asymptotically vanishes. Notably, its asymptotic behavior 
as well as limit laws for properly normalized number of particles at the source of branching, 
conditioned on the presence of particles at the origin, have different forms for d — 1, 2, 3, 4 and 
d > 5. Among the arising limit distributions one can find exponential and discrete ones along 
with a mixture of such laws. 

In the model of critical CBRW on Z d the behavior of number of particles located at an 
arbitrary point of the lattice remained unknown. The present work completes the picture. We 
study the asymptotic behavior in time of mean local particles numbers and that of probability 
of particles presence at a fixed point y ^ where = (0, 0, . . . , 0) G Z d . All the more, we 
obtain a conditional limit theorem for the properly normalized number of particles at such 
point y. It should be emphasized that in contrast to [4], [11], [13], [H] and [20] we admit 
the start of CBRW at an arbitrary point x G Z d and not only at the source of branching. 
Asymptotic properties of the number of particles at for CBRW with an arbitrary start point 
were investigated in [5J. 

The structure of the rest of the paper is the following. In section 2 we describe the model in 
detail and formulate three main results. Theorem [1] is proved in Section 3. Section 4 is devoted 
to construction of the auxiliary Bellman-Harris branching process. Thereupon we establish 
Theorems |2] and [3] in Section 5. 

2 Main results 

Now we dwell on the definition of a critical CBRW on Z d . At the initial time t = there 
is a single particle on the lattice located at a point x G Z d . If x ^ 0, the particle performs 
a continuous time random walk until the time of the first hitting the origin. The random 
walk outside the origin is symmetric, homogeneous, irreducible (i.e. a particle passes from an 
arbitrary u G Z d to any v G Z d with positive probability within a finite time) and has a finite 
variance of jumps. Accordingly, we assume this random walk be specified by an infinitesimal 
matrix A = (a(u,v)) UtV£Z d such that 

a(u,v) = a(v,u), a(u, v) = o(0, v — u) := a(v — u), u,v€l i d , 

Ea(v) = where a(0) < and a(v) > if v ^ 0, > ||-u|| 2 a(t;) < oo. 

If x = or the particle has just hit the origin it spends there an exponentially distributed 
time (with parameter 1). Afterwards, it either dies with probability a G (0, 1) producing before 
the death a random number of offsprings £ or leaves the source of branching with probability 
1 — a. In the latter case the intensity of transition from the origin to a point v ^ is given by 

a(0,v =-(l-aHfi(. 

0(0) 

At the origin the branching is determined by a probability generating function 

f(s):=^ = YTjks\ s €[0,1]. 

< • fc=0 



In [13] CBRW on Z is called critical if the following relations hold 

af(l) + (l-a)(l-h d ) = l and a 2 := a f" (1) < oo. (1) 

Here /i^ is the probability of the event that a particle leaving the origin will never return there. 
By the recurrence of a random walk on Z and Z 2 one has h\ = h 2 = 0. It is well known that 
h d G (0, 1) for d > 3. 

Newborn particles are located at the origin at the birth moment. They evolve according to 
the scheme described above independently of each other as well as of the parent particles. The 
number of particles located at a point y G Z d at time t > is denoted by //(£; y). 

The goal of the paper is three-fold. Firstly, we find the asymptotic behavior (as t — > oo) 
of the mean number of particles m(t;x,y) := E x /j,(t;y) located at a point y G Z d , y ^ 0, at 
time t > (everywhere the index x means that our CBRW starts at x G Z d ). Secondly, we 
retrieve the asymptotic behavior of the probability q(t; x,y) := P x (fj,(t; y) > 0) of the presence of 
particles at the point y at time t. Thirdly, we establish a limit theorem for properly normalized 
local numbers /i(£; y) conditioned on fi(t; y) > as t — > oo. 

To formulate the main results of the paper we introduce some more notation. Let p(t; x, y) 
be the transition probability from x to y within time t > for a random walk on Z d generated 
by matrix A. Set 

G x (x, y) := / e- xt p(t; x, y) dt, A > 0, x,y e Z d . 
Jo 

Note that the Green's function G (x,y) := lim A ^ + G\(x,y) is well-defined and takes finite 
values for d > 3 by virtue of the transience of our random walk on Z d , d > 3. One can check 
(see [133) that h d = (aGo(0,0)) _1 , d G N, where a := -a(0). 

As shown in [22], Theorem 2.1.1 (see also [13]), for any fixed x, y G Z d , one has 

P{t;x,y)~ ^, p(t;0,0) -p{t;x,y) ~ ^ 1+d/2 , t -^ oo, (2) 

where 7 d := ((27r) d |det ^(0)|)" 1/2 , 0(0) := £* 6 z*a(*> 0) cos(«, 0), G [-7r,7r] d , 

1 



to- ^ 



a^a^ 



7-W : = 7^TT- d I (v, zfeWo^/ 2 dv, z G Z d , 
o/ i,jG{i,...,4 lK - m > J ^- d 



and (•, •) stands for the scalar product in M. d . In particular, it follows that the value 

m d := 1 - (1 - a)a _1 + 2(1 - a)a' l G Q 2 {Q, 0) / tp(t; 0, 0) dt 

Jo 

is finite for d > 5. Set also g(s, i; x, j/) := 1 - ExS^* 5 ^, s G [0, 1], t > 0, x, y G Z d . For d = 2 we 
use the function 



J(s;y) := a 



POO 

/ (/(l-g(s,«;0,|/))-l + g(a,tt;0,2/))d«, s6[0,l], yeZ d . 
Jo 



The main results are contained in the following three theorems. For the sake of completeness 
their statements include the case y = studied earlier in [!]-[§], [H], [13], [14] and [20] . 



Theorem 1 Let x,y G Z . TTie following relations are valid being different for y ^ ami 
1/ = 0, namely, as t — > oo, 



7i 



m(t;x,y)~—p, m(t;x,0 



7i a 



vT 



m{t;x,y) 



72 



m(t; x, 



;i-a)V* 

72« 



G (a;,0)Go(0,y) 

"H*; a, y) t= — > m (^ a, °) 

G (z,0)G (0,y) , , 

m(r; x, y) ~ , m(t\ x, 0) 



(l-a)t 
aG (x,0)G (0,0) 



27773(1-0)^ 
aG (x,0)G (0,0) 



m{t;x,y) 



74 hit 
;i-a)G (a;,0)Go(0,y) 

aG 2 (0,0)m d 



m(t; x, 0) — t- 



74(1 — a) hit 
G (x,0) 



G (0,0)m d 



Theorem 2 For x,y E Z d and t — > 00 i/ie following formulae hold true 



q{t;x,y) 



2(1 -a) 
cr 2 7i a\/tlnt 



d=l 



gfoz^^h-^L-J^y)). // y 0. (/: 2 



g(t;x,0) ' 
q(t;x,y) 

q(t;x,y) ~ 



72 



1-J(0;0)), d = 2 
d = 3 



(1 -a)t 
4vr7 3 (l-a)G'o(x,0 



a 2 aG[](0,0)Vtlnt ' 
374(1 — a)Go(x, 0) hit 
a 2 aG[j(0,0)i ' 
2m d G (x,0) 



q(t;x,y) 



a 2 G (0,0)t 



d = 4 
d> 5 



d=l, 
d=2, 

d = 3, 

d = 4, 

d> 5. 



where for d = 2 and s G [0, 1] £/ie stnct inequalities J(s;y) < (1 — s)(l — a)/a, |/ ^ 0, and 
J(s; 0) < 1 — s are valid. 

Theorem 3 Given x, y G Z d , A G [0, 00) and s G [0, 1], one has, as t — > 00 ; 



limE x .^exp<| EMt . y) ^ t . y)>0) 



»(t;y) >o 



A + l 



d — 1, d = 3 or d > 5, 



t— >oo 



lim e. ( y"»» 1 m«; ») > 0) = (1 - a) ; : fy y)) . »*o, rf = 2 , 

Km E. ( «**> I „fe 0) > 0) = -'(JWOj-^iO)), d = 2 , 



i— >oo 



'-- L '' V " W I E a; (/x(t;y)Kt;j/)>0) 



1- J(0;0) 

12 2 



3 3 2 + 3A' 



d = 4. 



Observe that the normalizing factor E x (/j,(t;y)\/j,(t;y) > 0) arising in Theorem [3] is exactly 
m(t; x, y)/q(t; x, y) and the asymptotic behavior of the functions m{t\ x, y) and q(t\ x, y) is given 
by Theorems [1] and EJ, respectively. 



To establish Theorem [T] it is useful to invoke the forward and backward Kolmogorov's 
differential equations (considered in appropriate Banach spaces) for mean numbers of particles 
at different points of the lattice and also the resulting integral equations (see [22])- As for 
Theorems [2] and El note that for proving results in [I] , [TT] , [T3] , [H] and [20] concerning the 
number of particles at the origin the method of introduction of an auxiliary Bellman-Harris 
branching process with particles of two types was efficient. However, for proving Theorems [2] and 
[3] we have to involve a Bellman-Harris branching process with particles of six types. To apply 
the latter method we attend to a new notion of the hitting time with taboo in the framework 
of a (non-branching) random walk on Z d . More precisely, we use our recent results (see [8]) on 
the asymptotic behavior of the tail of the (improper) cumulative distribution function of this 
time. Due to that one can employ the theorems by V.A.Vatutin for Bellman-Harris branching 
processes with particles of several types (see, e.g., [16]-[T9]). Afterwards we have to deal with 
sophisticated analytic estimates of the solutions of the parametric integral equations (see, e.g., 
HE, H3], HI and TO- 

3 Proof of Theorem U 

Let us recall some useful results employed within this section. According to [ID] . Ch.3, Sec. 2, 
the transition probabilities p(t; x, y), t > 0, x, y G Z d , of the random walk generated by matrix 
A satisfy the backward Kolmogorov's equations 

(1 ' l[f:rlj) =(Ap(t;.,y))(x), p{0;x,y) = Sy{x). (3) 



dt 



d\ 



Here (Ap(t;-,y))(x) = ^2 zeZ dCi(x,z)p(t;z,y) and S y (-) is a column vector in the space Z 2 (^' 
with zero components except for the component 1 indexed by y. In a similar way, the backward 
Kolmogorov's equations for m(t;x,y), t > 0, x, y G Z d , (see, e.g., Theorem 2.1 in [23]) take the 
form 

dm ^ ,V) = (Am(t;;y))(x)+P c (A m(t;;y))(x), m{0;x,y) = 5 y (x), (4) 

where A = (a(u,v)) uv£Zd := A+(a~ 1 (l — a) — 1) A A, A := 5 5q (T stands for transposition) 
and (3 C := (1 — a)a~ 1 G^ 1 (0, 0). Here we follow the notation of 



Lemma 1 For each y G Z d , the function m(t; y, y) is non-increasing in t. 

Proof. The monotonicity of m(-;y,y) for SBRW on Z d was established in Lemma 3.3.5 of 
[22] . The key step of its proof was to use self-adjointness of the operator H := A + (3 C A where 
(3 C := G ( 5" 1 (0, 0). For CBRW the analog of H is the non self-adjoint operator H := A + /3 C A . 
However, Lemma 3.1 in [23] permits to pass to (self-adjoint) symmetrization of H and then 
apply Lemma 3.3.5 in [22J. Further argument is similar to the proof of Theorem in [2lJ. □ 

Equation (J3J) was obtained by differentiating at s = 1 the following backward Kolmogorov's 
equation for the generating function F(s,t;x,y) := E x s^ t ' y \ s G [0,1], t > 0, x, y G Z d , 
(see [23]) 

dF(s ^ y) = (AF(s, t; -, y)) (x) + (Aj(F(s, t; •, y))) (*), F(s, 0; x, y) = s^\ (5) 



Here f(s) := a(f(s) — s), s G [0, 1], is an infinitesimal generating function of the number of 
offsprings of a parent particle. We will employ fl5]) in Section 5. 

In Lemma [2] we derive a counterpart of the forward Kolmogorov's equation for the fuc- 
tion F(s,t;x,y), s G [0,1], t > 0, x, y G Z d , and, as a consequence, the forward Kol- 
mogorov's equation for m(t;x,y). Recall that A denotes an adjoint operator for A and 

\^A*m(t; x, •) J (y) = J2 z ez d m (*5 x ' z )^' 3/)- 

Lemma 2 For s G [0, 1], t > 0, x, y G Z d , £/ie following relation holds true 

9F{S ^ *' V) = {s-l) Y, a(z,y)E x s^ f ,(t ] z) + (s-l)a(y,y)E x s^- 1 ^t ] y) 

+ 6o{y)7(s)E x s^ t ^- 1 f i(t;y), F(s,0;x,y) = s 5 ^. (6) 

Moreover, one has 

dm{t d ^ y > = (Xm(t; x } •)) (y) + J c (A m(t; x, •)) (y), m(0; x, y) = 4(y)- (7) 

Proof. As usual in derivation of forward Kolmogorov's equations, we consider all possible evo- 
lutions of the particles population within the time interval [t, t + h) and let h — > 0+. To justify 
arising passages to the limit we involve the Lebesgue theorem on dominated convergence and 
useful estimates for transition probabilities (see proof of Lemma 3 in [ID], Ch.3, Sec. 2). We also 
benefit from finiteness of the mean total number of particles M(t;x) := E x (XLgz d M^i z )) f° r 
each x G Z d and t > 0. The latter observation is true since the last function belonging to /^(Z^) 
is a solution of the linear differential equation in fll]) with the initial condition M(0; x) = 1 for 
all x (instead of S y (x) in (01)), see [23] . 

Equation (JT]) is an immediate consequence of fl6j due to formula m(t; x, y) = d s F(s, t; x, y)| s =i- 
We also take into account that f'(l) = (3 C in view of ([I]). □ 

Consider equations (@| and (0) as inhomogeneous ones for differential equation ([3]) in Banach 
space l OD (7j d ). Applying the variation of constant formula (see [9], Ch.2, Seel) we infer that 

t 



m(t;x,y)=p(t;x,y) + (l / p(t — u; x, 0)m'(u; 0, y) du 

\ I- a J J 

+ — / p(t — u; x, 0)771(1/; 0, y) du, (8) 

1-a Jo 

m(t;x,y)—p(t;x,y) + I l) / m(t — u; x, 0)p'(u; 0, y) du 

+ /3 C m(t — u; x, 0)p(u; 0, y) du. (9) 

Jo 

An analogous result for SBRW on Z d can be found in [22], Theorem 1.4.1. Now we can give 
Proof of Theorem [Q To find the asymptotic behavior of m(t;x,y), t — >■ 00, x,y G Z d , 
y 7^ 0, we estimate each of the summands in the right-hand sides of (jH]) and when i^0 
and x = 0, respectively, as £ — > 00. Namely, we will show that, for d = 1 and d = 2, the main 
contribution to the asymptotic behavior of the right-hand side of ([SI), as well as of fl9J, is due 



to the first summand. However, for d > 3, the asymptotic behavior of the right-hand sides of 
(jHJ) and (Q is determined only by the third summands. It is worth mentioning that, for d = 1 
and d = 2, the third summands in (jSj) and (jUj) vanish in view of equality /3 C = 0. 

Let x = 0. The asymptotic behavior of the first summand in the right-hand side of (j5j) is 
given by fl2]). The estimate of the second summand could be obtained on account of Lemma 6 
in [13] and, in particular, relation (20). However, to avoid verifying the bounded variation of 
the functions p(t;x,y) and m(t;x,y) in variable t we choose another approach consisting in 
direct estimation of the second summand. Recall that representation (2.1.15) in [22] entails the 
inequalities p'(t; 0, 0) < 0, p'(t; 0, 0) < p'(t; 0, y), p"(t; 0, 0) > and p"(t; 0, 0) - p"(t; 0, y) > 0, 
t > 0. Then by virtue of (J2J) as well as the classical results on differentiating the asymptotic 
formulae (see, e.g., [3], Ch.7, Sec. 3), for d € N, one has 

//. n n N d ld ,, n , ,, v (d + 2)^ d (y) 

P W °> °) ~ -2#/2+T' p(t;0,0)-p(t;0,y) 2fd/2+2 , t -»- oo. 

Whence taking into account Lemma 5.1.2 in [22] ("lemma on convolutions") and the already 
proved assertion of Theorem [1] for x = and y = we deduce that, as t — Y oo, 

/ m(t-u;0,0)j/(u;0,y)du= / m{t -u; 0,0) (p'(u; 0, y) -p'(u; 0,0)) du 
Jo Jo 

+ / m(t- u; 0, 0)p'(u; 0, 0) rfw = m(t; 0, 0) - m(t; 0, 0) + o{m{t; 0, 0)) = o{m{t; 0, 0)). (10) 
Jo 

Combining relations (|2]), OU) and ( TTUj) we establish Theorem [TJ for d = 1 and d = 2 when x = 0. 
The statement of Theorem [1] for d > 3 and x = follows from formulae (|2D, © and (TTUj) by 
Lemma 5.1.2 in [22] and in view of Theorem [1] for the known case x = y = 0. 
Let x ^ 0. Similarly to the case x = we see that 

/ p(t — u; x, 0)m'(u; 0, y) du = / m(t — u; 0, y)p'(u;x, 0) du — o(m(t; 0, y)), t — >• oo. (11) 
Jo Jo 

Thus, the combination of (EJ), (JHI) and ( TTTi) proves Theorem [1] for d = 1 or d = 2 and i^0. 
For d > 3 and i^Owe estimate the third summand in (JSJ) with the help of Lemma 5.1.2 in 
[22] . relation ([2]) and the assertion of Theorem [1] for d > 3 and x = established above. D 

4 Auxiliary Bellman-Harris branching process 

Let us briefly describe a Bellman-Harris branching process with particles of six types. It is 
initiated by a single particle of type i = 1, . . . ,6. The parent particle has a random life-length 
with a cumulative distribution function (c.d.f.) Gi(t), t > 0. When dying the particle produces 
offsprings according to a generating function fi(s), s = (s\, . . . ,sq) G [0, l] 6 . The new particles 
of type j — 1, ... ,6 evolve independently with the life- length distribution Gj(t) and an offspring 
generating function fj(s). Let M := (d S] fi\g=(i ) ... : i)) i . =1 6 be the mean matrix of the process. 
The Bellman-Harris branching process is called critical indecomposable if the Perron root of M 
(i.e. eigenvalue having the maximal modulus) equals 1 and for some integer n all elements of 
M n are positive (see, e.g., [12], Ch.4, Sec. 6 and 7). Denote the number of particles of type j 



existing at time t by Zj(t), t > 0, j = 1, . . . , 6. Set Fj(t; s) = Ej ( ]^[ =1 s • J , z = 1, . . . , 6, 

£ > 0, s*G [0, l] 6 , where the index i means that the parent particle is of type i. In other words, 
Fi(t; s) is a generating function of the numbers of particles of all types existing at time t given 
that the process is initiated by a single particle of type %. 

Before demonstrating how an auxiliary Bellman-Harris process can be constructed in the 
framework of CBRW on Z d we have to introduce some notation. Recall that in [8] a new notion 
of a hitting time with taboo was proposed for a (non-branching) random walk on 7L d generated 
by matrix A. Namely, let t~ z , y,z G Z d , y ^ z, be the time spent by the particle (performing 
the random walk) after leaving the starting point until the first hitting y if particle's trajectory 
does not pass z. Otherwise (if particle's trajectory passes point z before the first hitting y), 
r~ z = oo. Denote by H~ yz (t), t > 0, the improper c.d.f. of r~ z given that the starting point of 
the random walk is x G Z d . 

Return to CBRW on Z d . In this section we assume that CBRW may start at the origin 
or at a fixed point y ^ 0. We divide the particles population existing at time t > into 
seven groups. The particles located at time t at the origin (respectively, at y) form the first 
(respectively, second) group having cardinality fi(t; 0) (respectively, fi(t;y)). Next consider at 
time t a family of particles labeled by a collection (u, v, w) of lattice points, its cardinality 
being n UtVtW (t). It consists of the particles which have left u at least once within time interval 
[0,£], upon the last leaving u have yet reached neither v nor w but eventually will hit v before 
possible hitting w. Our third group corresponds to (u,v,w) = (0,y, 0), the fourth to (y, 0,y), 
the fifth to (0, 0, y) and the sixth to (y, y, 0). The seventh group comprises the rest of particles 
not included into the above six groups. Note that the last group consists of the particles having 
infinite life-length since after time t they will not hit the origin any more. So, after time t these 
particles will not produce any offsprings and have no influence on the numbers of particles in 
other six groups. 

Now we can introduce an auxiliary Bellman-Harris process and use it for the study of CBRW 
on Z d . Consider a six-dimensional Bellman-Harris process having the following c.d.f. Gi and 
generating function fa, i — 1, . . . , 6, 

G 1 (t) = 1-e-*, f 1 (<?) = af(s 1 ) + (l-a)H- yjO (0)s 2 + (1 - o) (^,o(°°) -^ tfl o(0))*s 

+ (1 - a)H 00y (oo)s 5 + (1 - a)(l - H 0y0 (oo) - H 00y {oo)), 
G 2 (t) = l-e~ at , f 2 {s) = H- 0)y (0)si + (H- 0>y (o O )-H- 0ty (0))s A 

+ H y,yA°°) S S + i 1 ~ H V,0,y(°°) ~ H y,y,o(°°)), 

t^W - 77= — 7 — \ 7T 1 — 777' h\ s ) - s 2, <^A{t) - — — -, — : — — 7zt, U[s) — si, 

Ho^oi 00 ) ~ H o,y,o(°) H y,o,y(°°) ~ H y,o,y(°) 

ry /J.^ ^0,0,j/W x f^\ /~i lx\ V,V,0v') r /->\ 

° 5 ^ = u (r^v h(s) = s u G 6 (t)= r y ' , f 6 (s) = s 2 , 

H 0,0,y(°°) H y ,yA°°l 

where H~ (oo) := \\m t -^oo H~ yz {t). The symmetry and homogeneity of the random walk 
generated by matrix A imply identities Hq v0 = H~ 0y and Hq 0v = H~ , whence G3 = G4 
and G5 = Gq. It is not difficult to see that for the branching process constructed in this way 
one has (/i(£; 0), //(£; y), fi , y , (t), fx y ,o, y (t), Ho,o, y (t), fx y , y ,o(t)) = [Z l (t], ..., Z e (t)), t > 0. 

Observe that the introduced Bellman-Harris branching process with particles of six types is 
critical indecomposable. Indeed, it is an easy computation task to check that all entries of M 6 
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are positive. Furthermore, if Hq o (0) 7^ (that is a(0, y) > 0) then already all entries of M 4 
are positive. Hence, the constructed process is indecomposable. To verify its criticality note 
that in view of Theorem 3 in [S] one can rewrite the first relation in ([T]) as follows 

a/'(l) = 1 - (1 - a) (H ^ y (oc) + ^ {(X f ) \ ) ■ (12) 

Then by inspecting the explicit expression for the characteristic polynomial of the mean matrix 
M we deduce that it has the form 

det (M - kI) = k 2 (k - 1)R(k) 

where J is a unit matrix, k6C and 

R(k) := k 3 + k 2 {1 - «/'(l)) + k (1 - a/'(l) - (2 - a)H^ y (oo) - (1 - «)(^ o (0)) 2 ) 
+ (1 - a)(fl^ i0 (oo) - ^,o(0)) 2 - (1 - a)(H^ y (oo)) 2 . 

The polynomial R(k) has no real roots greater than 1 because R(l) > and R'(k) > for 
K > 1. In fact, due to identity (TT2j) we obtain the representation with strictly positive summands 

R(l) = (1 - if - 0>2/ (oo))((l - a)i/ 0A ,(oo) + 1) + (1 - a)F - y>0 (oo)(^ i0 (oo) - fl^ >0 (0)) 
(1 - a)g - i/[0 (cx))(g - yi o(o ) - ^ o - y[O (0)(l - g - , y (cx)))) | (1 - a)(g -, i0 (oo)) 2 

1 - Ho,0, y (°°) 1 - H 0,0,y{° o ) 

Moreover, if k > 1 then 

«'(«) = 3k 2 + 2/c(l - af'(l)) + 1 - a/'(l) - (2 - a)H- 0!y (oc) - (1 - «)(^, o (0)) 2 
> 3 - 2H 00 iy (oo) - ^ v , (0) > 0. 

Thus, the greatest positive real root of the characteristic polynomial of M is 1. Hence, by the 
Frobenius theorem (see, e.g., Theorem 2 in [T2], Ch.4, Sec. 5) 1 is the Perron root of M. So, 
the auxiliary Bellman-Harris process is critical. 

Denote by v = (t>i, . . . , v 6 ) and u = (ui, . . . , u 6 ) the left and right positive eigenvectors corre- 
sponding to the Perron root of M such that (u, 1 ) = 1 and (v,u) = 1 where l = (l,...,l)eM 6 . 
Taking into account ffT2|) we rewrite the components of u and v in the convenient form 



1 - H 0,0,y(°°) H 0,y,o(°°) Mo\ 

«1 = «4 = ^5 = Jj , U 2 = U 3 = U 6 = y — , (13) 

U U{\-OL)H^ yfi {oo) U(l-a)(H^ (oo)-H^ (Q)) 

V\ - — , V 2 — —J- — — , V 3 - — , [i.4: 

V V(l - H 00y {oo)) V 

U{\ - a)H 00 (00) 
v 4 = v 2 {H 0>y!0 {oo)-H 0y>0 (0)), v 5 = ^^ ' ^6 = v 2 H 00y (oo) (15) 

where the auxiliary variables U and V are defined by way of 

U : = 3(1 - H ^ y (oo) + #0^,0(00)), 

V:=3-2af\l)-(2- a )H ^ y (oo)^l-a)((H 0iyfi (oo)-H ^ o m^ 
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Using decomposition f(l — x) = 1 — f'(l)x + f"(l)x 2 /2 + o(x 2 ), x — > 0+, along with formulae 
(IT2|) - (IT51) and the definition of f(s) = (/i(s), • • • , fe(s)), it is not difficult to verify by standard 
calculations that 

x - (v, 1 - / f 1 - wxj J ~ 5x 2 , x -> 0+, where 5 := — i °y ■ (16) 



In the next two lemmas we apply theorems proved in papers [16]-|19] to the constructed 
six- dimensional Bellman-Harris branching process and then reformulate the obtained results 
for CBRW on Z d when d > 5. Common to these theorems are the conditions of criticality 
and indecomposability of the Bellman-Harris process which were established above. Another 
common condition on the behavior of the function x — (v,l — f(l — ux)) is fulfilled due to (TT6|) . 
However, various Vatutin's theorems involve different assumptions on the order of asymptotic 
decrease of the tails of Gk{-), k = 1, . . . , 6. It is worth to mention that such asymptotic behavior 
was established in [5], Theorem 3. Namely, our result for d < 5 corresponds to condition of 
Theorem 1 in [16j whereas the cases d = 6 and d > 7 meet the respective conditions of 
Theorem 3 in [19] and Theorem 2 in |18|. 



Lemma 3 Given i/gZ 5 , y ^ 0, for CBRW on Z 5 one has 

q{t; 0,y) = o (i~ 3 / 4 ) , q{t; y,y) = o (t^ 4 ) , t -)• oo. 

Proof. To apply Theorem 1 in [16] to the six-dimensional Bellman-Harris process constructed 
above for CBRW on Z 5 we verify the conditions of that theorem. According to the definition 
of G(-) = (Gi(-), . . . , Gq(-)) and by Theorem 3 in [S] for d = 5, the variable (5 in condition 2) of 
Theorem 1 in [TB] is equal to 3/2 whereas the function L\{t) in the same condition tends to a 
constant, as t — > oo. The validity of condition 3) of Theorem 1 in [16] is implied by Theorem 1 
in [T7] (for our process the function L\(t) in this theorem tends to 1/B, as t — > oo, in view of 
(1TB"]) ) combined with the definition of G(-) and Theorem 3 in [8] for d = 5. Thus, we may employ 
Theorem 1 in [16J. Taking into account Theorem 3 in [8] for d = 5 and formulae ( 1T3|) - (1T5|) 
we deduce from Theorem 1 in [16] that lim^oo Ej(s 2 \Z{t) ^ 0) = 1 for each S2 G [0, 1] and 
i = l,2 (as usual, Z(t) = (Z 1 {t),..., Z 6 (t)) and = (0, . . . , 0) G M 6 ). Setting s 2 = in the last 
relation one has Pi(Z 2 (t) > 0) = o(Pj(Z(t) 7^ )), as t — > 00. Moreover, examining the proof of 
Theorem 1 in [T6] we can show that for our Bellman-Harris process the slowly varying function 
L*(x) in the assertion of that theorem turns equivalent to l/y/B, asi-> 0+. Consequently, the 
indicated in [T6] formula (0.4) can be sharpened in our case, namely, the function Pj(Z(t) 7^ ) 
has an order of decreasing t _3//4 , as t — > 00. Whence by the connection between CBRW on Z 5 
and the auxiliary Bellman-Harris process we complete the proof. □ 

Lemma 4 In the framework of CBRW on Z d with d > 6 the following relations hold true for 

yeZ d ,y^0, 

2m d 2m d G (0,y) 

^;o, y )~— , g(*;y.y)~ ff2Go(0>0)t . *->~. 



Proof. Let us apply Theorem 3 in [19J to our Bellman-Harris branching process when d = 
6. To this end we verify whether all the conditions of Theorem 3 in [19] are satisfied. In 
view of ffTB"]) relation (6) in [19] is valid for our process and the function Li(n) in (6) tends 
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to l/B, as n — )■ oo. Equality (7) in [19] is also satisfied due to (6) in [19] in view of the 
definition of G(-) and Theorem 3 in [8] for d = 6. Now we may apply Theorem 3 in [19J. 
In particular, it follows that for each i — 1,2 the expressions lim^oo Pj(Zi(t) = 0|Z(t) 7^ 0) 
and lim^^oo Pj(Z 2 (t) = 0|Z(t) 7^ ) coincide, are positive and strictly less than 1. Consequently, 
P l (Z l (t) > 0) ~ P<(Z 2 (t) > 0), as t -> 00. The asymptotic behavior of g(i; 0, 0) = Pi(Zi(£) > 0) 
and q(t;y,0) = P 2 (Zi(£) > 0) can be found in [6], Lemmas 2 and 4. Thus, Lemma H] is proved 
for d = 6. 

For d > 7 we will employ Theorem 2 in [18]. Condition (6) of that theorem is valid due to 
Theorem 1 in [T7] (for our process, the function Li(t) in this theorem tends to l/B, as t — > 00) 
by virtue of the definition of G(-) and Theorem 3 in [8] for d > 7. The definition of Gk{-) and 
Theorem 3 in [8] for d > 7 also imply that L t dGk < 00 for each k — 1, . . . , 6. So, all the con- 
ditions of Theorem 2 in [18] are satisfied and it follows that lim^oo Pi(Zk(t) = 0\Z(t) 7^ ) = 
for each k = 1, . . . , 6 and i = 1, 2. Hence, P 4 (^i(t) > 0) ~ Pi(Z 2 (t) > 0), t ->■ 00. Notably, the 
asymptotic behavior of g(t; 0, 0) = Pi(Zi(t) > 0) and q(t;y, 0) = P 2 (^i(t) > 0) can be found 
in [6], Lemmas 2 and 4. Lemma H] is proved for d > 7. □ 

Concluding this section we derive an integral equation in function g(-; 0,y), y ^ 0, which 
is a counterpart of equation (2.6) in [TTJ for g(-;0, 0). Our integral equation will be essen- 
tially used for proving Theorem [2] when d = 4. Before formulating the corresponding state- 
ment we have to introduce some more notation. Let t z be the time spent by a particle per- 
forming a random walk generated by matrix A until the first hitting a point z G 7L d . In a 
similar way, r~ is the time spent by the particle after leaving the starting point of the ran- 
dom walk until the first hitting the point z. If the starting point of the random walk is z 
then the first hitting z means the first return to z. Denote by H X)Z (t) and H~ z (t), t > 0, 
the (improper) c.d.f. of t z and t~, respectively, given that the starting point of the ran- 
dom walk is x G 7L d . Obviously, H x>z (t) = G 2 * H~ z {t) for t > and x, z G Z d . Set also 
K(t) := a f'i^dit) + (1 - a)d * Hq (t) and h(s) := a{f(l-s) -l + f(l)s), s G [0,1]. Note 
that the function K(t) and the function K^(t), d G N, arising in [TTJ and [T3], coincide for each 
t > 0. Thus, Lemma 2.3 in [TT] and Lemma 11 in [T3] in which the asymptotic properties of 
c.d.f. Kd(t) and its density kd(t) are established, as t — > 00, may be applied to our function K. 

Lemma 5 For y G 7L d , j/^ 0, one has 

q(t- 0, y) = (1 - a)G 1 * (ff^,(t) - H , v (t)) + g(-; 0, y) * K(t) - %(•; 0, y)) * G^t). (17) 

Proof. Recall integral equations (see, e.g., [12], Ch.8, Seel) for probability generating func- 
tions F(t; s) := {Fi(t\ s ),..., F 6 (t; s)) of a six- dimensional Bellman-Harris process 

Fi(t;s) = Si(l-Gi(t))+ J fi(F(t-u;s)^ dG^u), t > 0, 8*6 [0,1], 2 = 1,..., 6. 

By setting here s = (1, 0, 1, 1, 1, 1) and substituting the explicit formulae for Gj, j = 3, 4, 5, 6, 
and fi, i = 1, . . . , 6, we get six integral equations in functions Fi(t) := Fj(t; (1, 0, 1, 1, 1, 1)), 
t > 0, i — 1, . . . , 6. Substituting the fourth and the sixth ones into the second equation and 
solving the obtained renewal equation in F 2 (-) we find 



F ^) = G 2*(i-H ^ (t)-H ^*J2Zo H ^y {t \ 

+ F 1 * G 2 * (fl^ >0 (0) + (iW) " ^,o(0)) * ^" =Q H* % >y (t)) 
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where H x ^ r (t) := G<i * H~ z r (t), t > 0, x, z, r G Z d , z ^ r. Now we substitute the last equation 
as well as the third and the fifth equations in functions Fi into the first one. After some 
algebraic transformations we obtain the following non-linear integral equation in function i*\ 

F x {t) = l-aG 1 *(l-f(F 1 (t)))-(l-a)G 1 *{l-F 1 (-))*Ho fi {t) 

(l-a)G 1 *(H- y (t)-H 0>y (t)) (18) 

provided that the following two equalities are valid 

E°° — k — — S. — v °° h 

k=0 H 0,y,0 * H y,y,0 * H y ^ y (t), H 0y (t) = H 0yfi * 2^ fc=Q #0,0,2/ W> 

for each t > 0. The first of them is true since any trajectory from to of a particle performing 
a random walk on 7L d either passes y exactly k times, k — 1, 2, . . . , or does not hit y until the 
first returning to 0. Similar argument justifies the second equality as well. Recall that due to 
the connection between the CBRW on Z d and the constructed Bellman-Harris process one has 
q(t;0,y) = Pi(Z 2 (t) > 0) = 1 — iq(t). Hence, rewriting ffTSl as an equation in q(t;0,y) we 
come to (fTTD. D 



5 Proofs of Theorems 2 and 3 



First of all, we derive some integral equations to be treated in this section. Consider equation ([5]) 
as inhomogeneous one for differential equation (j3]) in Banach space l 00 ('Z d ). By the variation of 
constant formula we infer (for a similar deduction see [5]) that 

q(s, t; x, y) = (1 — s)m(t; x, y) — \ m(t — u;x,0)h(q(s 1 u;0,y))du (19) 

Jo 

where q(s, t;x, y) — 1 — F(s, t;x,y), s G [0, 1], t > 0, x, y G 7L d . Substituting x = in the last 
equation we come to an integral equation in function q(s,t; 0,y) 



(20) 



q(s,t; 0, y) = (1 — s)m(t; 0, y) — \ m(t — u;0,0)h(q(s,u;0,y)) du. 

Jo 

Note that g(0, £; x, y) is equal to q(t; x, y). Thus, on account of (TT9l) one has 

q(t;x,y) = m(t;x,y)— / m(t — u;x,0)h(q(u;0,y)) du. (21) 

Jo 

Substituting x = in fT2T|) we derive an integral equation in function q(t; 0, y) 

q(t; 0, y) = m(t; 0, y) - [ m(t - u; 0, 0)h{q{u; 0, y)) du. (22) 

Now let us prove Theorems [2] and [3] for x = 0. Since their proofs depend ondGN essentially, 
we have to consider the cases d=l, d = 2, d = 3, d = 4 and d > 5 separately. Evidently, 
Theorem [2] for x = and d > 6 is implied by Lemma HI Due to Lemmas [TH3] and equation (T221 
the proof of Theorem [2] for x = in the respective cases d — 1, d — 2, d — 3 and d = 5 mainly 
follows the scheme proving, respectively, Theorem 2 in [20], Theorem 2 in [4], Theorem 4 in 
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[13] (item 3) and Theorem 4 in [13] (item 4). Moreover, by virtue of Lemma [5] the proof of 
Theorem [2] for x = and d — 4 is similar to that of Theorem 1.1 in [11]. So, we give only a 
few comments on the proof of Theorem [2] for x = and d < 5. 

If rf = 1 then the equality J °° h(q(u; 0, j/)) du = (1 — a)a _1 is valid. Furthermore, in view of 
©, ©, © and Theorem 5 in [13] one gets the useful estimate 

m(t; 0, y) - (1 - a)a _1 m(t; 0, 0) = O (t~ 3/2 ) , £ -> oo. 



a) a x . 



When d = 2 one can check the strict inequality J(0;y) = J °° h(q(u;0,y)) du < (1 
However, if d = 3 then J °° h(q(u; 0, y)) du = (1 — a)a~ 1 Go(0, y)GQ 1 (0, 0) and 

m(£; 0, y) - (1 - a)a -1 G (0, Z/)G^ 1 (0, 0) m(t; 0, 0) = 0(t _1 ), t -4- oo, 

by virtue of (jSJ), ([9]) and Theorem 5 in [13]. For d — 4 the first summand in (ITTj) is c^t -1 ), 
£ — ?• oo, by Lemma 3 in [6] and it does not contribute to the (main term of) asymptotic behavior 
of q(t; 0, y). As for d — 5, one has J °° h(q(u; 0, y)) dit = (1 — a)a _1 G (0, y)G ( 7 1 (0, 0) and 

m{t; 0, |/) - (1 - a)a -1 G (0, y)G l (0, 0) m(t; 0, 0) = O (r 3/2 ) , £ ^ oo, 

in view of (J2J), ([9]), Theorem 5 and Corollary 1 in [13J. Thus, Theorem [2] is proved for x = 0. 

Turn to Theorem [3] when x = 0. The proof of Theorem[3]for x = is similar to those of The- 
orem 4 in [T4], Theorem 2 in [4] and Theorem 4 in [H] for d — 1, 3, d — 2 and d > 5, respectively. 
Note only that the constant c* arising in the proof of Theorem[3]for x = in contrast to its coun- 
terpart in Theorem 4 in [H] is equal to a 2 7 2 a/(2(l-a)) and cr 2 aG 3 (0, 0)G (0, y)/(8vr 2 7 2 (l-a)) 
when d — 1 and <i = 3, respectively. At last, the constant c^ appearing in Theorem 4 in [6] 
equals (1 — a)Go(0,y)a 2 / (2aGo(0, 0)rri^) in the case of Theorem [3] for x = and d > 5. Since 
the limit theorem for /z(£; 0) when d = 4 was established by another approach, namely the 
moment method, we give the detailed proof of the limit theorem for /j,(t; y) when d = 4. So, to 
complete the proof of Theorem [3] for x = we dwell on the case d = 4 in detail. 

Set s(t) := s(t; A) = exp{-Aln 2 t/(c*t)} where c* := a 2 aG 3 (0, 0)G (0,y)/(37 2 (l-a)), t > 
and A > 0. By Theorems [1] and [2] for x = and d = 4 we see that 



m(t] 0, y) c*t 
g(t;0,y) ~lr7? 



Eo(Ai(*;y)lA«(*;y) > o) = T"' " - r^, *^°°- ( 2 3) 



The inequality 1 — e z < z for z > yields 



f { A In 2 £ w(w;y)l\ Aln 2 t Aln 2 t 

g(s(t), u; 0, j/) = E I 1 - exp <^ — V J < — ^— E /i(u; y) = — ^-m(u; 0, y) 

where u > and t > 0. By virtue of this estimate combined with Theorem 1 and the inequality 
h(z) < o" 2 2; 2 (being true for z > small enough) one has for t large enough 

nt/ In 3 t 

/ m(t — u; 0,0) h(q(s(t),u;0,y))du 

Jo 

a 2 A 2 ln 4 t /-' /ln3t „ N , N , Pl (t;A)lnt , s 

^ c , 2t2 y ™ 2 fa 0,y)m(t-w;0,0)d M = ^Y . (24) 
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Here p\ G U and U is the class of all bounded functions p(t; A) vanishing as t — > oo uniformly 
in A G [0, b], whatever positive b is taken. In a similar way, we obtain 

t 

m(t — u\ 0, 0)h(q(s(t), u; 0, y))du 

t-t/ In 2 t 

a 2 A 2 ln 4 £ /"* 2/ « n / ~ ~n , p 2 (t;A)lnt 
^ S7^ / m 2 (n;0,y)m(t-n;0,0)rfn=^^ (25) 

C l J t-t I In 2 1 t 

for p2 G W. It is not difficult to show that uniformly in u G [£/ In 3 t,t — t/ In 2 1] 

lnw~lnt, ln(t — w)~lnt, t — >• oo. (26) 

These facts, Theorem [TJ and the relation h(z) ~ a 2 z 2 /2, z — y 0, allow us to claim that 

^t-t/ln 2 t 
't/ln 3 t 

g 2 (s(t), M ;0,i/)ciM(l + p3(t;A)) 



/•t-t/ In 2 t 

J(t;A): = / h(q(s(t),u;0,y))m(t-w,Q,Q)du 

Jtl In 3 t 



a 2 m(t;0,0) r*-*J ln2 * 



2 



</ In 3 t 



where P3 G W. After changing the variable u = tv and using Theorems 1 and 2 for x = and 
(1 = 4 we get 

o rl-l/ln 2 t 

7 fo A ) = o ^ n W g 2 (s(t),tt;;0, 2 /)rfi;(l + p4(t;A)), ft eW. (27) 

2g(t; o,y) A/in 3 * 

In the last integral the function q(s(t; X),tv; 0,y) can be replaced by q(s(tv; Xv),tv; 0,y). In- 
deed, as 1 — e~ z < z for z > 0, we have 

|g(a(t; X),tv; 0, j/) - g(s(ti>; Xv),tv; 0, j/)| 

At; In 2 (it;) 1 / Aln 2 t 



= E (exp| — — jz(tu;y)j-exp| —p(tv;y) 

,. ( f A(-21ntlnt;-ln 2 t;) . N ]\ A(-21ntlnt; - ln 2 i;) , 
< E (l-exp| ^ — J -pt(tv;y)jj <^ — J -m(tv;0,y). 

Since functions z\nz and zln z are bounded for z G (0, 1), by virtue of Theorems [1] and [2] for 
x = along with relation (1261) we see that uniformly in v G [1/ In t, 1 — 1/ In t] and < A < A 
with an arbitrary positive A 

q(s(t;X),tv;0,y) q(s(tv;Xv),tv;0,y) 

7 n 7 n * U, t — ¥ OO. (28) 

q(tv;0,y) q(tv;0,y) 

Set <p(t;A) := g(s(t; A),*; 0,y)/(Xq(t; 0,j/)), t > 0, A > 0. Then dividing both sides of flU by 

Xq(t; 0, |/) and using (T2~4"]) -( |2"8]) along with Theorem [2] for x = and relation 1 — e~ z ~ 2, z — > 0, 

we obtain 

3A /' 1 - 1 /in 2 t 
<p(t; A) = 1 + p 5 (£; A) - — / (f 2 (tv;Xv)dv, p 5 eU. 

* Jl/\n 3 t 
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Changing the variable w = Xv leads to the following relation 

3 ^(l-l/la 2 ') „ ftw 
'A/ In 3 t 

The argument similar to the proof of Theorem 4 in [H] establishes that 



y?(£;A) = l + p 5 (£;A) -- / y> T i» aw. 

2 Jx/inH \ X / 



lim<p(£;A)=<p(A) = — ^— , 0<A<A , (29) 

t-xx> SX + 2, 

where Ao is some positive number and <p(X) is the unique solution of the equation 

3 f x 
tp(\) = 1 - - / tp 2 (w) dw, A > 0. 
2 Jo 

Invoking the definition of cp(t\ A) we rewrite relation fl29|) by way of 



lim E < exp 



A In 2 £u(£;y)1 1 4 , , , s 1 2 



t-K» I ^ C*t 

for < A < Ao- Since both the Laplace transform of a non-negative random variable and the 
function 1/3 + 2/3 • 2/(3A + 2) are analytic and bounded in the domain {A : Re X > 0} C C, 
by the uniqueness theorem for analytic functions relation ( 130]) is valid for each A with Re A > 
(for an analogous deduction see, e.g., [IS]). Combining ( 123]) and (130]) we complete the proof of 
Theorem [3] for x = and d = 4. Thus, Theorem [3] is proved for x = 0. 

Next we prove Theorems [2] and [3] when i / 0. As a preliminary we derive some more 
integral equations. In the framework of CBRW on Z d , the parent particle can either hit the 
point or not within time interval [0,t\. In the latter case at time t there is a single particle 
on Z d located at the point y or outside it. Consequently, 

E x s^ y) = E x s^I(t <t) + E xS ^I(t > t, n(t; y) = 1) + E x s^I(t > t, /x(t; y) = 0) 
= E x s^ y) I{r o <t) + sP x {T O >t,fi{t;y) = l) + P x {T O >t,fi{t;y) = 0) (31) 

where !(•) stands for the indicator of a set. Evidently, the first summand in (13T]) can be rewritten 
in the form 

E x s^I(r < t) = J s^dP x = J E x (s^\r )dP x 

{r <t} {r <t} 

= [ E x (s^\t = u) dH Xi0 (u) = f E s^- U ^ dH X;0 (u). (32) 
Jo Jo 

It is easily seen that the probability at the second summand in (J3"T1) can be represented as 
follows 

fc=0 ^to*(l-G 2 (£)) when x + y, (33) 
P y (r >t,^t;y) = l) = ^=0^0 * i 1 " G 2 (t)). (34) 
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It also turns convenient to write the third summand in fl3T|) in the form 

k=Q H* y k yfi *(l-G 2 (t)) if x^y, (35) 
P y (r o >t^(t;y) = 0) = 1 - H y , (t) - YT k=Q Ky,o * (1 - G*(t)). (36) 

Combining relations f l3~Tj) -( l3"6"]) we come to the desired integral equations 

q(s, t; x, y) = (1 - s)^,^ * $^ fc=0 ^*° * ^ ~ ^^ + / g ( S ' t ~ M; °' y ) dH ^ o( ^> if x ^' 
q(s, t- y, y) = {l- s) J2Z H *v k yfi * t 1 " ^(*)) + / <K«, * - «; 0, y) dH yfi (u). 
In particular, for s = one has 

g(t; a;, y) = #^ j0 * ^" =q #** i0 * (1 - G 2 (t)) + f q{t - w; 0, y) dfl a , i0 («) when x ^ y, (37) 
q(t;y,y) = y^ H* y k y0 *(l-G 2 (t))+ / g (t- U ;0,j,)^,o(«). (38) 

Now we have the tools for proving Theorems [2] and [3] for x ^ 0. To establish Theorem [2] 
for x 7^ and d ^ 2 we employ equations (|3"7|) and (1381) . It is not difficult to see that the first 
summands in the right side of (|37|) and (|38|) are equal to pit; x, y) — f Q p(t — u; 0, y) dH X)0 {u) 
for x 7^ y and x = y, respectively. The latter expression can be rewritten as follows 

p(t;x,y)- / p(t-u;0,y)dH xfi (u) 
Jo 

= p{t; x, y) - p{t; x,0)+ f {p{t - u; 0, 0) - p{t - u; 0, y))dH x p(u) (39) 

Jo 

due to the obvious relation p(t; x, 0) = f Q pit — u; 0, 0) dH Xt0 (u). The asymptotic behavior of 
the first summand at the right-hand side of (I39p is given by formula ([2]) whereas the asymptotic 
behavior of the second summand in (1391) can be found with the help of relation (J2J) , Lemma 3 
in [6] and Lemma 5.1.2 in [22] . Finally, the first summands in fl37|) and (138]) are 0(t~ 3 ' 2 ) when 
d—1 and 0{t~ d / 2 ) when d > 3. Hence, the first summands in (j3"7|) and (1381) are o(g(t; 0, y)), 
as £ — > oo, by Theorem [2] for x = 0. Moreover, on account of Lemma 3 in [6] and Lemma 5.1.2 
in [22] we reveal that the last summands in (13T1) and (1351) are equivalent to g(t; 0, y) and 
g(t; 0,|/) G (x, 0)G " 1 (0, 0) for d = 1 and d > 3, respectively, as t — > oo. Hence Theorem [2] is 
proved for x ^ and d ^ 2. As for Theorem [2] when i^0 and d = 2 as well as Theorem [3] 
for x 7^ 0, we only note that their proofs bear on analysis of equations ( 1191) and (I2~TI) . Since 
the proofs are similar to that of Theorem 5 in [6], they are omitted. So, Theorems [2] and [3] are 
proved completely. □ 

The author is grateful to Associate Professor E.B.Yarovaya for permanent attention and to 
Professor V.A.Vatutin for valuable remarks. Special thanks are to Professors I.Kourkova and 
G. Pages for invitation to LPMA UPMC (Paris-6) where this work was started. 
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